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We study the potential energy between static charges in G2 gluodynamics in three and four dimensions. Our 
work is based on an efficient local hybrid Monte-Carlo algorithm and a multi-level Luscher-Weisz algorithm 
with exponential error reduction to accurately measure expectation values of Wilson- and Polyakov loops. Both 
in three and four dimensions we show that at intermediate scales the string tensions for charges in various G2- 
representations scale with the second order Casimir In three dimensions Casimir scaling is confirmed within 
one percent for charges in representations of dimensions 7, 14, 27, 64, 77, 77', 182 and 189 and in 4 dimensions 
within 5 percent for charges in representions of dimensions 7, 14, 27 and 64. In three dimensions we detect 
string breaking for charges in the two fundamental representations. The scale for string breaking agrees very 
well with the mass of the created pair of glue-lumps. 



There is compelling experimental evidence that the fundamental constituents of QCD, quarks and gluons, never show up as 
asymptotic states of strong interaction - rather they are confined in mesons and baryons. Understanding the dynamics of this 
confinement mechanism is one of the challenging problems in strongly coupled gauge theories. There are convincing analytical 
and numerical arguments to believe that confinement is a property of pure gauge theories (gluodynamics) alone and that the 
underlying mechanism should not depend on the number N of colours. Confinement is lost at high temperatures and for gauge 
groups with a non-trivial center the trace of the Polyakov loop 



vanishes in the confined low-temperature phase and is close to an element of the center in the deconfined high-temperature 
phase. In gluodynamics or gauge theories with matter in the adjoint representation the action and measure are both invariant 
under center transformations, whereas the Polyakov loop transforms non-trivially and hence serves as order parameter for the 
global center symmetry. This means that the center symmetry is realized in the confined phase and spontaneously broken in the 
deconfined phase. 

In the vicinity of the transition point the dynamics of the Polyakov loop is successfully described by effective Sd scalar field 
models for the characters of the Polyakov loop [1^]. If one further projects the scalar fields onto the center of the gauge group 
then one arrives at generalized Potts models describing the effective Polyakov-loop dynamics [5]. The temperature dependent 
couplings constants of these effective theories have been calculated ab initio by inverse Monte Carlo methods in [3]. 

With dynamical quarks in the fundamental representation the center symmetry is explicitly broken and the Polyakov loop 
points always in the direction of a particular center element. In a strict sense the Polyakov loop ceases to be an order parameter. 
This is attributed to breaking of the string connecting a static 'quark anti-quark pair' when one tries to separate the charges. It 
breaks via the spontaneous creation of dynamical quark anti-quark pairs which in turn screen the individual static charges. 

The pivotal role of the center for confinement also follows from a recent observation relating the Polyakov loop with center 
averaged spectral sums of the Dirac operator [6-8]. More precisely, for gauge groups with non-trivial center one can relate the 
expectation value of the Polyakov loop to dual condensates. This result could finally explain why for gauge groups with a non- 
trivial center and fundamental matter the transition temperatures for the deconfinement and chiral phase transitions coincide. On 
the contrary, for gauge theories with adjoint matter the two transition temperatures can be very different [9, 10]. 

To clarify the relevance of the center for confinement it suggests itself to study pure gauge theories whose gauge groups have 
a trivial center. For such theories the string connecting external charges can break via the spontaneous creation of dynamical 
'gluons' such that the Polyakov loop acquires a non-vanishing expectation value for all temperatures, similarly as it does in QCD 
with dynamical fermions. Here the simple gauge group 50(3) suggests itself and indeed the 5*0(3) gauge theory has been 
studied in great detail on the lattice, see for example [11]. Unfortunately, via the non-trivial first homotopy group tti (50(3)) = 
Z2 the lattice gauge theory 'detects' its simply connected universal covering group SU{2). To avoid the resulting lattice artifacts 
one should investigate theories with simply connected gauge groups with trivial center. 
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TABLE I. Centers Z of simple lie groups. 



group 
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Dr, r odd 
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center Z 
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From Tab. I, taken from [12], one reads off that the smallest simple Lie group with these properties is the 14-dimensional 
exceptional Lie group G2. This is one reason why the group in Bern investigated G2 gauge theories with and without Higgs 
fields in series of papers [13-15]. In their pioneering works it has been convincingly demonstrated that G2-gluodynamics shows 
a first order finite temperature phase transition without order parameter from a confining to a deconfining phase. In this context 
confinement refers to confinement at intermediate scales, where a Casimir scaling of string tensions has been reported [16]. On 
large scales strings will finally break due to spontaneous gluon production and the static inter-quark potential is expected to 
flatten [17]. However, the threshold energy for string breaking in G'2-gauge theory is rather high and all previous attempts to 
detect this flattening have been without success. In the present paper we shall demonstrate that string breaking for charges in 
the fundamental and adjoint representations of G2 takes place at the expected scales. To that aim we implemented a shghtly 
modified Ltischer-Weisz multistep algorithm for high-precision measurements of the static inter-quark potential. 

The present paper deals with G'2-gluodynamics in 3 and 4 dimensions. The simulations are performed with an efficient and 
fast implementation of a local HMC algorithm. Below we shall calculate the potentials at intermediates scales for static charges 
in the 7, 14, 27, 64, 77, 77', 182 and 189-dimensional representations. We show that in 3 and 4 dimensions the string tensions on 
intermediate scales are proportional to the second order Casimir of the representations. The high-precision measurements in 3 
dimensions confirm Casimir scaling within 1 percent. In 4 dimensions Casimir scaling for the lowest 4 representations is fulfilled 
within 5 percent. In 3 dimensions we also calculated the static potential for widely separated charges in the two fundamental 
representations. In both cases we see a flattening of the potential which signals the breaking of the cormecting string. The energy 
where string breaking sets in is in fuU agreement with the independently calculated masses of the glue lumps formed after string 
breaking. 

n. THE GROUP G2 

The exceptional Lie-Group G2 is the automorphism group of the octonion algebra or, equivalently, the subgroup of SO{7) that 
preserves any vector in its 8-dimensional real spinor representation. This means that the 8-dimensional real spinor representation 
of Spin(7) branches into the trivial representation and the 7-dimensional fundamental representation of G2- The 14-dimensional 
fundamental representation of G2, which at the same time is the adjoint representation, arises in the branching of the adjoint 
of SO{7) according to 21 — )• 7 ® 14. The 27-dimensional representations of SO{7) acting on symmetric traceless 2-tensors 
remains irreducible under G2. In this work we need the foUowing branchings of S'0(7)-representations to G2: 

7^7, 21^^14 8 7, 27-)-27, 35^-27 7 1, 77^^77. (2) 

For explicit calculations it is advantageous to view the elements of the 7-dimensional representation of G2 as matrices in the 
defining representation of SO{7), subject to seven independent cubic constraints [15]: 

Tabc = Tdef Qda deb Qfc- (3) 

Here T is a total antisymmetric tensor given by 

Tl27 = 2^154 = T'leS = ^235 = ^264 = ^374 = T576 = 1. (4) 

The gauge group SU{2i) of strong interaction is a subgroup of G2 and the corresponding coset space is a sphere [18], 

G2lSU{2,) ~ S^. (5) 

This means that every element U of G2 can be factorized as 

U^S-V with V e ^[/(S) and 5 e G2/5'?7(3), (6) 

and we shall use this decomposition in our simulations. The short exact sequence 

= 774 (S*^) ^ 7r3(5J7(3)) ^ 7r3(G2) ^ 773(5*^) = (7) 

shows that 773(02) = Z and hence there should exist (j2-instantons of any integer topological charge. In the charge fc-sector 
there are at least 3fc magnetically charged defects [12]. 
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TABLE II. Representations of G2 with corresponding dimension and Casimir values. 



representation Tl 


[1,0] 


[0,1] 


[2,0] 


[1,1] 


[3,0] 


[0,2] 


[4,0] 


[2,1] 


dimension d-jz 


7 


14 


27 


64 


77 


77' 


182 


189 


Casimir eigenvalue C-r, 


12 


24 


28 


42 


48 


60 


72 


64 


Casimir ratio 


1 


2 


7/3 


3.5 


4 


5 


6 


16/3 



Any irreducible representation of G2 is characterized by its highest weight vector /x which is a linear combination of the 
fundamental weights, = + 9M(2)> with non-negative integer coefficients p, q called Dynkin labels. The dimension of an 
arbitrary irreducible representation 7?. = [p, g] can be calculated with the help of Weyl's dimension formula and is given by 

dn = dimp,g = ^(1 + + g)(2 + p + (7)(3 + p + 2<7)(4 + p + 3g)(5 + 2p + 3g). (8) 

Below we also use the physics-convention and denote a representation by its dimension. For example, the fundamental repre- 
sentations are [1,0] = 7 and [0, 1] = 14. However, this notation is ambiguous, since there exist different representations with 
the same dimension. For example [3, 0] — 77 and [0, 2] = 77' have the same dimension. An irreducible representation of G2 
can also be characterized by the values of the two Casimir operators of degree 2 and 6. Below we shall need the values of the 
quadratic Casimir in a representation [p, q], given by 

Cn = Cp.q = 2p2 + 6g2 + 6pg + lOp + 18g. (9) 

For an easy comparison we normaUze these 'raw' Casimir values with respect to the defining representation by ^ = Cp^q/Cifi. 
The normaUzed Casimir values for the eight non-trivial representations with smallest dimensions are given in Tab. 11. 

Quarks and gluons in G2 are in the fundamental representions 7 and 14, respectively. To better understand G2-gluodynamics 
we recall the decomposition of tensor products of these representations, 

7(8)7=107® 14 ©27 
7 (8 14 = 7 ® 27 ® 64 
14® 14= 1© 14 ©27 ©77 ©77' (10) 
70707=l©4-7®2-14®3-27®2-64®77' 
14 O 14 O 14 = 1 ® 7 ® 5 • 14 ® 3 • 27 ® • • • 

The decompositions (10) show that, similarly as in QCD, two or three quarks or two or three gluons can buUd colour singlets - 
mesons, baryons or glueballs. Since three gluons can screen the charge of a single (static) quark, 

7 14 14(8)14= 1®--- , (11) 

one expects that the string between two static quarks will break for large charge separations. The two renmants are two glue- 
lumps - charges screened by (at least) 3 gluons. The same happens for charges in the adjoint representation. Each adjoint charge 
can be screened by one gluon. 



Construction of characters from tensor products 



The character xn = tr 7^ of any irreducible representation 7?. is a polynomial of the characters X7 and Xi4 of the two fundamental 
representations 7 and 14. For example, the first two decompositions in (10) imply 

X27 = X7 • X7 - Xl - X7 - Xl4 

X64 = X7 • Xl4 - X7 - X27 = X7X14 - X7 + Xl + Xl4 

and yield the characters of the representations 27 and 64 as polynomials of X7 ^ind xu- From further tensor products of 
irreducible representions one can calculate the polynomial in xn =Po1k(x7, X14) for ai^Y irreducible representation TZ. For a 
fast implementation of our algorithms we also need reducible representations. In particular we use 

(7(8 7)s, (7(8 7(8 7)s, (7 (8 7 (8 7 (8 7)s, (7(8 7)^(8 14 (13) 
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where the subscript 's' denotes the symmetrized part of the respective tensor product. Comparing the reduction of representations 
for SO{7) and G2 and mapping representations from SO{7) to G2 the following characters of reducible representations can be 
computed 

X(7®7)s = X27 + Xli 
X(7(8>7(8i7), = X77 + X7, 

(14) 

X(7®7(8)7(8i7), = Xl82 + X77 + X27 + X64 + 2 Xl4 + X7, 
X(7(gi7),(g)14 = X189 + X27 + Xl- 



ni. CASIMIR SCALING AND STRING BREAKING FOR SU{N) GAUGE THEORIES 

In QCD quarks and anti-quarks can only be screened by particles with non- vanishing 3-aUty, especially not by gluons. Thus, in 
zero-temperature gluodynamics the potential energy for two static color charges is linearly rising up to arbitrary large separations 
of the charges. The potentials for charges in a representation TZ can be extracted from the 2-point correlator of Polyakov loops 
or the expectation values of Wilson loops with time-extent T according to 

(P7z(O)PK(i?0> = , (VK^(i?,T)) =e««-^^-(«). (15) 

With dynamical quarks the string should break at a characteristic length r;, due to the spontaneous creation of quark anti-quark 
pairs from the energy stored in the flux tube connecting the static charges. However, for intermediate separations r < n the 
string cannot break since there is not enough energy stored in the flux tube. 

For pure gauge theories we expect the following qualitative behavior of the static potential: At short distances perturbation 
theory applies and the interaction is dominated by gluon exchange giving rise to a Coulomb-Uke potential, V ^ —a/r, the 
strength a being proportional to the value C-jz of the quadratic Casimir operator in the given representation TZ of the charges; at 
intermediate distances, from the onset of confinement to the onset of color screening at r^, the potential is expected to be linearly 
rising, V ~ ar, and the corresponding string tension is again proportional to the quadratic Casimir; at asymptotic distance scales 
(partial) screening sets in such that the string tension typically decreases and only depends on the A^^-ahty of the representation. 
In particular for center-blind color charges or gauge groups without center the potential flattens. The characteristic length n 
where the intermediate confinement regime turns into the asymptotic screening regime is determined by the masses of the 
debris left after string breaking. The Casimir scaling hypothesis, according to which the string tension at intermediate scales 
is proportional to the quadratic Casimir of the representation [19], is exact for two dimensional continuum and lattice gauge 
theories and dimensional reduction arguments support that it also holds in higher dimensions. Within the Hamiltonian approach 
to Yang-Mills theories in 2 + 1 dimensions the following prediction for the string tensions has been derived [20] 

.94 

(^n = -;— CuC-R.. (16) 

For pure SU(2) and SU(3) gauge theories in three and four dimensions there is now conclusive numerical evidence for Casimir 
scaling from Monte-Carlo simulations: for SU (2) in 3 dimensions [19, 21] and in 4 dimensions [22-25] as well as for SU{3) in 
4 dimensions at finite temperature [26] and zero temperature [27-30]. In particular the simulations for SU{3) gluodynamics in 
[29] confirm Casimir scaling within 5% for separations up to 1 fm of static charges in representations with Casimirs (normalized 
by the Casimir of {3}) up to 7. String breaking for charges in the adjoint representation has been found in several simulations: 
In 3-dimensional SU (2)-gluodynamics with improved action and different operators in [31, 32] and in 4-dimensional SU (2)- 
gluodynamics in [33] with the help of a variational approach involving string and glueball operators. For a critical discussion of 
the various approaches we refer to [34], where string breaking in a simple setting but with an improved version of the Liischer- 
Weisz algorithm has been analyzed and compared with less sophisticated approaches. There is a number of works in which a 
violation of Casimir scaling on intermediate scales have been reported. For example, it has been claimed that in 4-dimensional 
S'i7(A^)-gluodynamics with larger A/" = 4, 6 the numerical data favor the sin-formula, as suggested by supersymmetry, in place of 
the Casimir scaUng formula [35]. The differences between the Casimir scaling law and sin-formula are tiny and it is very difficult 
to discriminate between the two predictions in numerical simulations. Indeed, in [36] agreement with Casimir scaling and sin- 
formula in 4-dimensions and disagreement in 3-dimensions has been claimed. In addition the high precision simulation based on 
the Liischer-Weisz algorithm in [37] point to a violation of the Casimir scaling law in 3-dimensional SU (2) gluodynamics. In a 
very recent paper Pepe and Wiese [38] reanalyzed the static potential for SU (2)-gluodynamics in 3 dimensions with the help of 
the Liischer-Weisz algorithm and confirmed Casimir scahng at intermediate scales and 2-ality scaling at asymptotic scales. 

For gauge theories with matter we expect a similar quaUtative behavior: a Coulomb-like potential at short distances, Casimir 
scaling at intermediate distances and (partial) screening at asymptotic distances. The string tension at asymptotic scales depends 
both on the iV-alities of the static color charges and of the dynamical matter. In particular, if dynamical quarks or scalars can 
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form center blind composites with the static charges then the potential is expected to flatten at large separations. To see any 
kind of screening between fundamental charges requires a full QCD simulation with sea quarks, which is demanding. Thus 
the earlier works dealt with gauge theories with scalars in the fundamental representation. For example, in [39] clear numerical 
evidence for string breaking in the 3-dimensional 5[/(2)-Yang-Mills-Higgs model via a mixing analysis of string and two-meson 
operators has been presented. Probably the first observation of hadronic string breaking in simulation of QCD3 with two flavors 
of dynamical staggered fermions using only Wilson loops have been reported in [40, 41]. Despite extensive searches for colour 
screening in 4-dimensional gauge theories with dynamical fermions the results are still preliminary at best. First indications for 
string breaking in two-flavor QCD, albeit only at temperatures close to or above the critical deconfinement temperature, have 
been reported in [42]. More recently Bali et al. used sophisticated methods (e.g. optimized smearing, improved action, stochastic 
estimator techniques, hopping parameter acceleration) to resolve string breaking in 2-flavor QCD at a value of the lattice spacing 
a^^ « 2.37 GeV and of the sea quark mass slightly below [43]. By extrapolation they estimate that in real QCD with light 
quarks the string breaking should happen at r?, « 1.13 fm. 

To measure the static potential and study string breaking three approaches have been used: correlations of Polyakov loops 
at finite temperature, variational ansaetze using two types of operators (for the string-like states and for the broken string state) 
and Wilson loops. Most results on Casimir scaling and string breaking have been obtained with the first two methods. This is 
attributed to the small overlap of the Wilson loops with the broken-string state. To measure Polyakov or Wilson loop correlators 
for charges in higher representations or to see screening at asymptotic scales one is dealing with extremely small signals down 
to 10~^". In order to measure such small signals one needs to improve existing algorithms considerably or/and use improved 
versions of the Liischer-Weisz multistep algorithm. 

For gauge groups with trivial centers like G2 , F4 or Eg, the flux tube between static charges in any representation will always 
break due to gluon production. The potential flattens for large separations and expectation values of the Polyakov loop never 
vanish [13]. However, for G2 it changes rapidly at the phase transition temperature and is very small in the low-temperature 
confining phase, see Fig. 1. Similarly as in QCD we characterize confinement as the absence of free colour charges in the 
physical spectrum [16, 44]. 
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FIG. 1. Phase transition on a 16 x 6 lattice in terms of the Polyakov loop in the fundamental representation. 



IV. ALGORITHMIC CONSIDERATIONS 
A. Local hybrid Monte-Carlo 

In simulations of gauge field theories different algorithms are in use. For 5[/(iV)-gluodynamics heat-bath algorithms based on 
the Cabibbo-Marinari SU{2) subgroup updates, often improved by over-relaxation steps, have proven to be fast and reliable. 
For QCD with dynamical fermions a hybrid Monte-Carlo (HMC) scheme is preferable. Based on [45] also local versions of 
HMC algorithms are available where single links are evolved in a HMC style. According to [46] the cost for the local hybrid 
Monte-Carlo (LHMC) is about three times more than for a combined heat-bath and overrelaxation (HOR) scheme for the case 
of 5[/(A^)-gluodynamics. 

For the exceptional gauge group G2 there exists a modification of the heat-bath update [13] which combines the heat-bath 
update for a 5[/(3)-subgroup with randomly distributed G2 gauge transformations to rotate the SU{3) subgroup through G2- 
In the present work we instead use a LHMC algorithm for several good reasons: First, the formulation is given entirely in terms 
of Lie-group and Lie-algebra elements and there is no need to back-project onto 6*2. The autocorrelation time can be controlled 
(in certain ranges) by the integration time in the molecular dynamics part of the HMC algorithm. Furthermore, one can use 
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a real representation of G2 and relatively simple analytical expressions for the two involved exponential maps to obtain a fast 
implementation of the algorithm. Finally, the inclusion of a (normalized) Higgs field is straightforward and does not suffer from 
a low Metropolis acceptance rate (even for large hopping parameters). 

The LHMC algorithm has been essential for obtaining the results in the present work. Since we developed the first implemen- 
tation for G2 it is useful to explain the technical details for this exceptional group. As any (L)HMC algorithm for gauge theories 
it is based on a fictitious dynamics for the link-variables on the gauge group manifold. The "free evolution" on a semisimple 
group is the Riemannian geodesic motion with respect to the Cartan-Killing metric 

ds% ^ KlT {dUU^^ ® dUU^^) . (17) 

In the fictitious dynamics the interaction term is given by the Yang-Mills action of the underlying lattice gauge theory and hence 
it suggests itself to derive the dynamics from the Lagrangian 

^ = ^ (^^.^Z^.'i)' - ^ymM, (18) 

where 'dot' denotes the derivative with respect to the fictitious time parameter r and 

SyuP] = ^ E (2 A^c - Ux,,u - Ul^,) (19) 

^ X.flU 

is the Wilson action. The Lie algebra valued fictitious conjugated link momentum is given by 

dL dL ■ 

Vx,^. = i^. -^=iU,,^,-p— = -iUx,^U-l, (20) 

and via a Legendre transform yields the pseudo-Hamiltonian 

H=^Y.^T^l^ + SYMm. (21) 

X,fX 

The equations of motion for the momenta are obtained by varying the Hamiltonian. The variation of the Wilson action 5ym[^] 
with respect to a fixed link variable Ux.fi is given by the corresponding staple variable Rx.fj., the sum of triple products of 
elementary link variables closing to a plaquette with the chosen link variable. Hence we obtain 



(22) 



The variational principle implies that the projection of the term between curly brackets onto the Lie algebra 02 vanishes, 

^X,IJ, — Pi^,x\g^- (23) 

Choosing a trace-orthonormal basis {Ta} of 02 the equations for the (L)HMC dynamics can be written as follows, 

'^x,^. = J2^r{FxMTa and U^,^ = ^^x,^Mx,^, (24) 

a 

with the "force" F^^fj. defined in (22). Now a LHMC sweep consists of the following steps: 

1 . Gaussian draw of the momentum variable on a given Unk. 

2. Integration of the equations of motion for the given link. 

3. Metropolis accept/reject step, 

4. Repeat these steps for all links of the lattice. 

This local version of the HMC does not suffer from an extensive SH oc V problem such that already a second order symplectic 
(leap frog) integrator allows for sufficiently large timesteps St. In condensed form the integration for a hnk variable yields 

U{t + 6t) = cxp {i^{t + 5T/2)5T)U{t). (25) 

For a large range of Wilson couplings /? in our simulations an integration length of T = 0.75 with a step size of 6t = 0.25 is 
optimal for minimal autocorrelation times and a small number of thermalisation sweeps. Acceptance rates of more than 99% are 
reached. Nevertheless, the most time consuming part of the calculations involves the exponential maps. A calculation for G2 
can be implemented fast and exact up to a given order in 5t as will be shown in the next section. 
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B. The exponential map ~^ 

For an efficient and fast computation of the exponential map we use the real embedding of the SU (3)-representation 3 ® 3 into 
G2, given by 



/l \ 



V(W) = fl^ 



Q€G2, with W€SU{3). 



W 
\0 W*J 

One can choose the unitary matrix to have block diagonal form with fin = 1. A possible choice for fl is 



(26) 



'1 ' 

,0 VQ, 



with Q = 



/O 1\ 

1 

1 
1 
1 

\0 1 0/ 



V 



V2 \i 1^ 



(27) 



Every element of G2 can be factorized as 



U^S-V(W) with SeG2/SU{3). 



(28) 



For a given timestep 5t in the molecular dynamics this factorization will be expressed in terms of the Lie algebra elements with 
the help of the exponential maps, 



exp{(5Tu} = expj^Ts} • expl^ro} with generators u e 02, OS V*(su(3)) 
fulfilling the commutation relations 

[o,D']=t)", [0,5] =5' and [s,s']=X)' +s". 



(29) 



(30) 



The generators s are orthogonal to the generators of the really embedded S'[/(3)-subgroup. To simplify the notation we absorb 
the time step 6t in the Lie algebra elements. 

The last exponential map in (29) can be calculated with the help of the embedding (26) and the exponential map for SU{3), 
W = cxp(rD), which follows from the Cayley-Hamilton theorem for S'C/(3)-generators, see [47]. The result can be expressed 
in terms of the imaginary eigenvalues wi, W2, ^3 of It) and the differences 5i = W2 — W3, 82 = W3 — wi and S3 = wi — W2 as 
follows: 



W = exp(rD) 



6i525i 



(ail + a„tt) + a„2it)^) 



(31) 



with expansion coefficients 



(32) 



i=l 



wherein one identifies w^+i and Wi. 
For the generators {Ui, . . . , U14} of G2 we use the real representation given in [44]. The su(3)-subalgebra formed by the 

elements {ui, .... Us} generates the really embedded 3 3 of SU (3) and the remaining generators {Ug, . . . , U14} generate the 
coset-elements S in the factorization (28). With this choice for the generators the real embedding (26) reads 



V(W) = 



'1 

.0 v^^ 



/ CI33 —^33 ^32 —^32 — &31 a3l\ 

^33 O33 ^32 (I32 031 631 

023 — &23 0'22 — &22 — &21 0-21 

623 O23 ^^22 Ct22 0,21 ^21 

&13 ai3 &12 ai2 ail &11 

\ai3 -613 ai2 -bi2 -hii an) 



(33) 
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where the entries are the real and imaginary parts of the elements of the S'J7(3)-matrix, Wij = aij + ibij. 
Finally, to parametrize the elements of the coset space we calculate the remaining exponential map 

S = exp {s} with s = ^ SjUs+i . (34) 

The result depends on the real parameter a = \\s\\ and the 6-dimensional unit- vector s = s/||s ||. In a 1 x 6-block notation the 
map takes the form 

^ _ I cos 2a — sin 2a 
I sin 2a s S± 

with 6-dimensional matrix 

S± = cos cr 1 + sin cr s_L + (cos 2(7 — COS a) ss^ + (1 — cos a)vv~^ . (36) 

The matrix s± is the 6x6 right-lower block of s in (34). The unit-vector = (s2, — si, S4, —§3, —sq, §5) defining the last 
projector in (36) is orthogonal to the unit-vector s defining the projector ss^. 

In the numerical integration we need the exponential map for elements u in 02- They are related to the generators used in the 
factorization by the Baker-CampbeU-Hausdoff formula, 

5tu = St{s + vi) + ^St'^[s,\3]-\ (37) 

Depending on the order of the symplectic integrator we must solve this relation for 5 and t) up to the corresponding order in St. 
For a second order integrator used in this work this can be done analytically since the commutator [s, o] does not contain any 
contribution of the sub-algebra su(3). The integrator used in the (L)HMC algorithm must be time reversible. It can be checked 
that time reversibility holds to every order in this expansion. To summarize, for a second order integrator the approximation (37) 
may be used in the exponentiations needed to calculate V and S. This approximation leads to a violation of energy conservation 
which is of the same order as the violation one finds with a second order integrator. In comparison to the exponentiation via 
the spectral decomposition the method based on the factorization (28) is more than ten times faster. It is also much faster than 
computing the exponential map for S0{7) via the Cayley-Hamilton theorem. 



C. Exponential error reduction for Wilson loops 

In the confining phase the rectangular Wilson loop scales as W{L, T) cx cxp(— crL -T). In order to estimate the string tension 
a we probe areas LT ranging from up to 100 and thus W will vary by approximately 40 orders of magnitude. A brute force 
approach where statistical errors for the expectation value of Wilson or Polyakov loops decrease with the inverse square root of 
the number of statistically independent configurations by just increasing the number of generated configurations will miserably 
fail. Nevertheless, convincing results on Casimir scaling on intermediate scales for representions with relative Casimirs 
< 5 have been obtained in [16] with a variant of the smearing procedure. When reproducing these results we observed 
that the calculated string tensions depend sensitively on the smearing parameter^ Thus to obtain accurate and reliable numbers 
for the static potential and to detect string breaking we implemented the multi-step Liischer-Weisz algorithm with exponential 
error-reduction for the time transporters of the Wilson-loops [48]. With this method the absolute errors of Wilson lines decrease 




This is not the case for the ratios of string tensions. 
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exponentially with the temporal extent T of the line. 
This is achieved by subdividing the lattice into rit 
sublattices l^i , . . . , containing the Wilson loop 
and separated by time slices plus the remaining sub- 
lattice, denoted by V, see figure on the right. At the 
first level in a two-level algorithm the time extent of 
each sublattice Vn is 4 such that rit is the smallest 
natural number with Ant > T+2. In the figure on the 
right T — 14 and the lattice is split into four sublat- 
tices Vi, V2, V3, V4 containing the Wilson loop plus 
the complement V. The Wilson loop is the product 
of parallel transporters W = 7^' 7^' 71 737271. If a 
sublattice V,i contains only one connected piece of 
the Wilson loop (as Vi and V4 do) then one needs to 
calculate the sublattice expectation value 



(7ri)?i 



1 



VU Tn e 



(38) 
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if Vn contains two connected pieces (as Vi and V3) 

then one needs to calculate (71i ® T^n- The updates in each sublattice are done with fixed link variables on the time-slices 
bounding the sublattice. Calculating the expectation value of the full Wilson loop reduces to averaging over the hnks in the 
nt + 1 time slices. 



{w) = (c((ri)i(r2 ® T^2 ■ ■ ■ {Tn,-i ® r^_i)„,-i(r„,>„,)) 



boundaries 



(39) 



(W) 



2nt-2 

tr n c^-: 



boundaries 



(40) 



Here C is that particular contraction of indices that leads to the trace of the product W — T2 ■ ■ ■ Tn^-i%it%n-i • • • 7171. In a 
two-level algorithm each sublattice Vn is further divided into two sublattices Vn.i and Vlt.2, see right panel in the above figure, 
and the sublattice updates are done on the small sublattices Vn,k with fixed link variables on the time slices separating the 
sublattices Vn,k- This way one finds two levels of nested averages. Iterating this procedure gives the multilevel algorithm. 

Since the dimensions dn grow rapidly with the Dynkin labels [p, q\ - for example, be- 
low we shall verify Casimir scaling for charges in the 189-dimensional representation 
[2,1]- it is difficult to store the many expectation values of tensor products of parallel 
transporters. Thus we implemented a slight modification of the Liischer-Weisz algo- 
rithm where the lattice is further split by a space slice with hyperplane orthogonal to 
the plane defined by the Wilson loop, see figure on the right. The sublattice updates 
are done with fixed link variables on the same time slices as before and in addition 
on the newly introduced space slice. Instead of nt sublattices containing the Wilson 
loop we now have 2nt — 2 sublattices. But now every sublattice contains only one 
connected part of the Wilson loop and (39) is replaced 
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VI 
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VI3 
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— I— 
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Vli 
— I— 



V12 



V 
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V^i2 



Til 



V02 



V21 



'13 
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An iteration of this procedure by additional splittings of the time slices leads again 
to a multilevel algorithm. In the present work we use a two level algorithm with 
time slices of length 4 on the first and length 2 on the second level. We calculate 
{W) for Wilson loops (and hence transporters 71) of varying sizes and in different 
representations. To avoid the storage of tensor products of large representations we 
implemented the modified algorithm as explained above. 

We also applied the Liischer-Weisz algorithm to calculate the correlators of two Polyakov loops {P-ji(Q)Pti{R)) on larger 
lattices. In this case the complete lattice is divided into sublattices separated by time slices, hence there is no complement V. 
Since the Polyakov loops are only used for lower-dimensional representations we have not split the lattice by a spatial slicing 
but used tensor products similar to Eq. (39). Actually for the calculations of Polyakov loop correlators we used the three-step 
Liischer-Weisz algorithm. 
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V. STRING TENSION AND CASIMIR SCALING IN G2 GLUODYNAMICS 



The static inter-quark potential is linearly rising on intermediate distances and the corresponding string tension will depend on 
the representation of the static charges. We expect to find Casimir scaling where the string tensions for different representations 
TZ and TV scale according to 

^ = ^ (41) 

with quadratic Casimir c-jz- Although all string tensions will vanish at asymptotic scales it is still possible to check for Casimir 
scaling at intermediate scales where the linearity of the inter-quark potential is nearly fulfilled. 

To extract the static quark anti-quark potential two different methods are available. The first makes use of the behavior of 
rectangular Wilson loops in representation TZ for large T, 

{Wn{R,T))=ei^p{nn{R)-VniR)T) with Vn{R) = m - + (t-rR- (42) 

The potential can be extracted from the ratio of two Wilson loops with different time-extent according to 

We calculated the expectation values of Wilson loops with the two-level Ltischer-Weisz algorithm and fitted the right hand 
side of (43) with the potential Vtz{R) in (42). The fitting has been done for external charges separated by one lattice unit up 
to separations R with acceptable signal to noise ratios. From the fits we extracted the constants 77^, a-ji and cr-ji entering the 
static potential. For an easier comparison of the numerical results on lattices of different size and for different values of (3 we 
subtracted the constant contribution to the potentials and plotted 

Vn{R) = Vn{R) - m (44) 
in the figures. The statistical errors are determined with the Jackknife method. In addition we determined the local string tension 

^ I P\ _ Vn{R + p)-Vn{R) 

moc,TZ \R+ 2) p ' ^ 

given by the Creutz ratio 

^ I PA 1 {Wn{R + P,T)){Wn{R,T + T)) _ an , ^ 

a^o.,n {Wn{R + p,T + r)){Wn{R,T)) " R(rTp} + ^'^^ 

The second method to calculate the string tensions uses correlators of two Polyakov loops, 

Vn{R) = In {Pn{0)Pn{R)) ■ (47) 
Pt 

The correlators are calculated with the three-level Liischer-Weisz algorithm and are fitted with the static potential Vji{R) with 
fit parameters jn, ctn and g-r.- Now the local string tension takes the form 

1 , {PTz{0)Pn{R + p)) 



A. Casimir scaling in 3 dimensions 

Most LHMC simulations are performed on a 28^ lattice with Wilson loops of time-extent T = 12. To extract the static potentials 
from the ratio of Wilson loops in (43) we chose r = 2. The fits to the static potential (42) for charges in the fundamental 
7 -representation and for values /3 = 30, 35 and 40 yield the lattice parameters a, 7 and a given in Tab. III. To check for scaUng 
we plotted the potentials in 'physical' units, V/ p, with mass scale set by the string tension in the 7-representation, 

/X = V^, (49) 

as function of pR in Fig. 2. We observe that the potentials for the three values of /? are the same within error bars. In addition 
they agree with the potential (in physical units) extracted from the Polyakov loop on a much larger 48"^ lattice. 



11 



TABLE III. Potential for charges in the 7-representation. 





/3 = 30, L = 28 


/3 = 35, L = 28 


/3 = 40, L = 28 


^ = 30, L = 48 




0.185(8) 


0.160(4) 


0.147(5) 


0.197(1) 


a 


0.0881(7) 


0.0752(3) 


0.071(4) 


0.098(1) 




0.046(1) 


0.0340(8) 


0.024(1) 


0.0435(3) 




FIG. 2. Continuum scaling of the fundamental potential. 



The fitted constants a-jzj^n and a-jz of the potential (42) for the eight smallest representations are given in Tab. IV. The 
Casimir scaling of coefficients becomes apparent when they are divided by the corresponding coefficients of the static potential 
in the 7-representation. 

The local string tensions extracted from the Creutz ratio can be determined much more accurately as the global string tensions 
extracted from fits to the static potentials. Tab. V contains the local string tensions for static charges in the eight smallest 
representations for p — I and different R in (46), divided by the corresponding local string tensions in the 7-representation. The 
results are insensitive the the distance R in the Creutz ratio. They agree within 1 percent with the values for the Casimir ratios 
C!ji — C-ji/C^ given in the last row of that table. 

In Fig. 3 we plotted the values for the eight potentials V7, . . . , Visg (with statistical errors) measured in 'physical units' /i 
defined in (49). The distance of the charges is measured in the same system of units. The linear rise at intermediate scales is 
clearly visible, even for charges in the 189-dimensionaI representation. 

Fig. 4 contains the same data points rescaled with the quadratic Casimirs of the corresponding representations. The eight 
rescaled potentials fall on top of each other within error bars. This implies that the full potentials for short and intermediate 
separations of the static charges show Casimir scaling. 

To further check for Casimir scaling we calculated the local string tensions with p = 1, this time for all R between 1 and 
10 and not only for i? = 0, 1, 2 as in Tab. V. The horizontal lines are the values predicted by the Casimir scaling hypothesis. 
Clearly we see no sign of Casimir scaling violation on a 28^-Iattice near the continuum at /3 = 40. Of course, for widely 



TABLE IV. Fit-parameters of static potentials. 



7^ 


7 


14 


27 


64 


77 


77' 


182 


189 




0.147(5) 


0.29(1) 


0.34(1) 


0.51(1) 


0.58(1) 


0.74(2) 


0.83(1) 


0.77(2) 




0.147 


0.145 


0.146 


0.146 


0.145 


0.148 


0.138 


0.144 


mln 


1 


1.97 


2.31 


3.46 


3.94 


5.03 


5.64 


5.23 


an 


0.071(4) 


0.145(8) 


0.16(1) 


0.24(1) 


0.27(1) 


0.36(1) 


0.37(1) 


0.36(1) 




0.071 


0.0725 


0.069 


0.069 


0.068 


0.072 


0.062 


0.068 




1 


2.04 


2.25 


3.38 


3.80 


5.07 


5.21 


5.07 




0.024(1) 


0.048(2) 


0.057(3) 


0.086(4) 


0.099(5) 


0.120(6) 


0.157(6) 


0.132(6) 




0.024 


0.024 


0.024 


0.025 


0.025 


0.024 


0.026 


0.025 




1 


2.00 


2.37 


3.58 


4.12 


5.00 


6.54 


5.50 
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TABLE V. Scaled local string tension. 



7^ 


7 


14 


27 


64 


77 


77' 


182 


189 


a,j(l/2)/<77(l/2) 


1 


1.9996(3) 


2.3327(5) 


3.4981 


3.997(2) 


4.996(3) 


5.991(5) 


5.328(4) 


aTC(3/2)/<77(3/2) 


1 


1.99897 


2.3311 


3.495(5) 


3.994(4) 


4.9897 


5.991 


5.321(9) 


^^(5/2)7(77(5/2) 


1 


1.9961 


2.3271 


3.484(5) 


3.9807 


4.961 


5.94(2) 


5.291 




1 


2.0000 


2.3333 


3.5000 


4.0000 


5.0000 


6.0000 


5.333 
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FIG. 3. Unsealed potential with /? = 40 on a 28^ lattice. 



separated charges in higher dimensional representations the error bars are not negligible even for an algorithm with exponential 
error reduction. 



B. Liischer term 



In Tab. IV we have seen that the dimensionless coefficient a-ji in the static potential scales with the quadratic Casimir, similarly 
to the string tension. The corresponding term, if measured at distances where the flux tube has akeady formed, is referred 
to as Liischer term. Its value has been calculated by Liischer for charges in the fundamental representation, in d dimensions 
a = {d — 2)7r/24, and it is believed to be universal [49]. The value a — 7r/24 in 3 dimensions is off the results in Tab. III. 
However, since the coefficients in this table are fitted to the static potential from i? 1 to values of R with acceptable signal to 
noise ratio, they contain contributions from the short range Coulombic tail. To calculate a-ji at intermediate distances we better 
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FIG. 4. Scaled potential with /? = 40 on a 28^ lattice. 
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FIG. 5. Ratio of the local string Tension with /? = 40 scaled on a 28 lattice for the eight smallest representations. 



use the (local) Liischer term 



i?' In iPn{0)Pn{R + P)) {Pn{0)Pn{R - p)) 



{PniO)Pn{R)) {PTziO)PTz{R)) 



(50) 



with p = 1. In Fig. 6 we plotted the local Liischer term for charges in the 7-representation on a larger 48'^-lattice with [3 = 30. 
Our data at intermediate distances are in agreement with the theoretical prediction aj = 7r/24 « 0.131. 




FIG. 6. Local Liischer term on a 48^ lattice at /? = 30. 



C. String breaking and glue-lumps in 3 dimensions 



To observe the breaking of strings connecting static charges at intermediate scales when one further increases the separation of 
the charges we performed high statistics LHMC simulations on a 48"^ lattice with (3 — 30. We calculated expectation values of 
Wilson loops and products of Polyakov loops for charges in the two fundamental representations of G2 ■ When a string breaks 
then each static charge in the representation TZ at the end of the string is screened by N{TZ) gluons to form a colour blind glue 
lump. We expect that the dominant decay channel for an over-stretched string is string — > gluelump + gluelump. For a string 
to decay the energy stored in the string must be sufficient to produce two glue-lumps. According to (11) it requires at least 3 
gluons to screen a static charge in the 7-representation, one gluon to screen a charge in the 14-representation and two gluons to 
screen a charge in the 27-representation. We shall calculate the separations of the charges where string breaking sets in and the 
masses of the produced glue-lumps. The mass of such a quark-gluon bound state can be obtained from the correlation function 



'N{n) 



CAT) 



n=l 



N(Tl) 
n=l 



cx exp(-TO7^r), 



(51) 



14 




where TZiUyx) is the temporal parallel transporter in the representation TZ from x to y of length T. It represents the static 
sources in the representation TZ . The vertical line means projection of the tensor product onto that linear subspace on which the 
irreducible representation TZ acts, 

(14(»14(»•••(g)14) = 7^®••• . (52) 
For example, for charges in the 14-representation the projection is simply 

F^,{x) =F^^{x), where F^^T^ ^ F^,. (53) 

14, a 

For charges in the 7-representation we must project the reducible representation 14 (g) 14 (X) 14 onto the irreducible representation 
7. Using the embedding of G2 into SO{7) representations one shows that this projection can be done with the help of the totally 
antisymmetric e-tensor with 7 indices, 

F^.{x) F^,{x) ® F^,{x) cx FP,{x)F^^,{x)F;;,{x)eabcdefgTPjlTJg. (54) 

7, a 

Fig. 7 shows the logarithm of the glue-lump correlator (51) as function of the separation of the two lumps for static charges in 
the fundamental representations 7 and 14. The linear fits to the data yield the glue-lump masses 

TO7 = 0.46(4), mi4 = 0.767(5). (55) 

Thus we expect that the subtracted static potentials approach the asymptotic values 

Vji — !■ 2rnK - jn- (56) 

With the fit-values 77 = 0.197(1) and 714 = 0.381(2) we find 

VV/m — ^3.46 , Vu/ii — ^ 5.52. (57) 

Fig. 8 shows the rescaled potentials for charges in the fundamental representations together with the asymptotic values (57) 
extracted from the glue-lump correlators. Within error bars both potentials flatten exactly at separations of the charges where the 
energy stored in the flux tube is twice the glue-lump energy. 

A good approximation for the string breaking distance is then given by Vji{R") ~ 2 m-ji. Assuming Casimir scaling for the 
coefficients a-jz , -f-n and a-ji in the static potential we obtain 

4(?-"»rK7-*))' 

Inserting the result from the last row in Tab. Ill and the glue-lump masses we find /ii?" — 4.00 and fiRf^^ = 3.28. These values 
agree well with the separations fiR in Fig. 8 where the static potentials flatten such that string breaking sets in at scales predicted 
by formula (58). Fig. 9 shows the local string tensions in the two fundamental representation and Fig. 10 their ratios. Especially 
the last plot makes clear that the string connecting charges in the adjoint representation break earlier than the string connecting 
charges in the 7-representation. The formula (58) predicts Rf^ = 9.40 and just above this separation the ratio of local string 
tensions cri4(i?) / cr-j{R) shows indeed a pronounced knee. 
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FIG. 8. Potential for both fundamental representations (lattice size 48 , /3 — 30) and corresponding glue-lump 




FIG. 9. Local string tension (48^ lattice, /3 = 30). 




FIG. 10. Casimir scaling of local string tension (48'' lattice, /3 = 30). 
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TABLE VI. Parameters of the quark anti-quark potential in 4 dimensions. 





/3 = 9.7, L = 14 


^ = 10, L = 14 


/3 = 9.7, L = 20 




0.83(8) 


0.74(4) 


0.68(9) 


a 


0.40(7) 


0.33(3) 


0.28(8) 




0.07(2) 


0.042(9) 


0.11(1) 



D. Casimir scaling in 4 dimensions 

In this last section we present our results for the static potential in 4 dimensions. The local HMC-simulations have been 
performed on a small 14^ and a larger 20'' lattice for different values of (5. The static potentials and local string tensions have 
been extracted from (43) and (46), where the expectation values have been calculated with a two-step Liischer-Weisz algorithm. 
Tab. VI contains the fits to the parameters in the potential for static charges in the 7-representation for these lattices and values 
for/3. 

Fig. 1 1 shows the static potentials in 'physical units' ^ = y/o^ for charges in the 7, 14, 27 and 64-dimensional representations 
and coupling /3 = 9.7 as function of the distance between the charges in physical units. The corresponding value for aj is taken 
from Tab. VI. The same coupling has been used in [16] on an asymmetric 14"^ x 28 lattice. After normahzing the potential with 
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FIG. 11. Unsealed potential at /? = 9.7 on a 14 lattice. 

the quadratic Casimirs they are identical within error bars, as can be seen in Fig. 12. Our findings are in complete agreement 
with the results in [16] on Casimir scaling in 4-dimensional G'2-gluodynamics at /? = 9.7 and our accurate results on Casimir 
scaling on intermediate scales in 3-dimensional G2-gluodynamics. 
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FIG. 12. Scaled potential at /3 = 9.7 on a 14"' lattice. 



Figs. 13 and 14 show the corresponding results for a weaker coupling /3 = 10 closer to the continuum limit. For this small 
coupling we can measure the potential only up to separations /ii? w 1.5 of the charges. But we can do this with high precision 
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and for higher-dimensional representations. As for /3 = 9.7 we find that the potentials normalized with the second order Casimirs 
fall on top of each other. This confirms Casimir scaling for G'2-gluodynamics in 4 dimensions for charges in representations 
with dimensions 7, 14, 27, 64, 77, 77', 182 and 189. 
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FIG. 13. Unsealed potential at /3 = 10 on a 14" lattice. 
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FIG. 14. Scaled potential at /? = 10 on a 14* lattice. 



Finally we simulated on a much larger 20'^ lattice at (3 — 9.7 in order to calculate the static potential for larger separations of 
the static quarks. Unfortunately the distance /ii? « 3 is still not sufficient to detect string breaking, see Fig. 15. But again the 
potentials normalized with the quadratic Casimirs shown in Fig. 16 are equal within error bars. 
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FIG. 15. Unsealed potential at /? = 9.7 on a 20^* lattice. 
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FIG. 16. Scaled potential at ^ = 9.7 on a 20"* lattice. 



TABLE VII. Fit-parameters of static potentials (20'' lattice, /3 = 9.7). 



TZ 7 14 27 

jTza 0.68(9) 1.39(4) 1.61(3) 
ina/C'T^ 0.68 0.695 0.690 

an 0.28(8) 0.60(2) 0.69(2) 
am/C'n 0.28 0.30 0.295 

aW 0.11(1) 0.21(1) 0.251(9) 
(Jna^/Cn 0.11 0.105 0.107 




In Tab. VII we have listed the fit-values for the parameters of the potentials on the larger 20'* lattice for static charges in 
the representations with dimensions 7, 14 and 27. For all representation we find Casimir scaling of all three parameters in the 
potential. Unfortunately the fit-parameters cannot be determined reliably in the 64-representation with the present data. This 
is attributed to larger errors for the potentials at intermediate scales, see Fig. 15, so that the parameters can only be determined 
from the ultraviolet part of the potential for this representation (R < 3) which is rather Coulomb-Iike than linearly rising. Much 
more conclusive are the local string tensions calculated on the larger lattice (now up to the 64-representation). Tab. VIII contains 
the local string tensions divided by the local string tensions in the 7-representation. These normalised values are constant up to 
separations of the charges where the statistical errors are under control. Compared to the corresponding numbers in 3 dimensions, 
see Tab. V, we now see a slight dependence of the local string tensions from Eq. 45 on the distance R. Despite of the lower 
precision of the results in 4 dimensions compared to the corresponding results in 3 dimensions we again confirm Casimir scaling 
on intermediate scales within 5 percent. 

All our simulation results for the local string tensions an{R) normalized by ar{R) on a U^-lattice with /? e {9.7, 10} and 
on a 20'*-lattice with /3 = 9.7 and for /ii? < 1.5 are collected in Fig. 17. The horizontal lines in this figure show the prediction 
of the Casimir scaling hypothesis. The normalized data points are compatible with each other and with the hypothesis. 



TABLE VIII. Scaled local string tension (20* lattice, /3 = 9.7). 



TZ 7 14 27 64 

(7^(1/2)7(77(1/2) 1 1.973(1) 2.294(1) 3.396(8) 

(7,j(3/2)/(77(3/2) 1 1.987(3) 2.303(4) 3.44(2) 

(7k(5/2)/(77(5/2) 1 1.92(1) 2.28(3) — 

Cn 1 2.0000 2.3333 3.5000 
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FIG. 17. Scaled local string Tension with P = 9.7, 10 on 14'', 20* lattices. 
VI. CONCLUSIONS 

In the present work we implemented an efficient and fast LHMC algorithm to simulate G2 gauge theory in three and four 
dimensions. With only a slight modification we can include a (normalized) Higgs field in the 7-representation. The corresponding 
results for the phase diagram of G'2-Yang-Mills-Higgs theory will soon be presented in a companion paper The algorithm has 
been optimized with the help of the coset decomposition of group elements and the analytic expressions for the exponential 
maps for the two factors. In addition we implemented a slightly modified Liischer-Weisz multi-step algorithm with exponential 
error reduction to measure the static potentials for charges in various G2-representations. The accurate results in 3 dimensions 
show that all parameters of the fitted static potentials show Casimir scaling, see Tab. III. The global string tensions extracted 
from these fits show that possible deviations from Casimir scaling, if they exist, must be less than 4 percent. We also extracted 
the local string tensions from the Creutz ratios to obtain even more precise data. This way we confirm Casimir scaling with 1 
percent accuracy. Thus we conclude that in 3-dimensional G2-Gluodynamics the string tensions show Casimir scaling for all 
charges in the representations with dimensions 7, 14, 27, 64, 77, 77', 182 and 189. In passing we can check the scaling formula 
(16) for the string tension cr-jiiP) as function of the coupling /3 cx 1/g^ [20]. On a fixed lattice this formula implies that the 
product fPaTi{l3) should be independent of /3. Using the values for the string tension aj in Tab. Ill we obtain 





;3 = 30 


/3 = 35 


/3 = 40 


(77 


0.046(1) 


0.0340(8) 


0.024(1) 




41.4 


41.7 


38.4 



The numbers in the last row show that the scaling cr oc l//3^ is almost fulfilled. In the present work we did not attempt to further 
clarify this interesting point by simulating at many /3-values and using the more accurate local string tensions. 

For charges in the two fundamental representations we performed LHMC simulations on larger lattices to detect string break- 
ing at asymptotic scales. In 3 dimensions we observe that string breaking indeed sets in at the expected scale where the energy 
stored in the flux tube is sufficient to create two glue lumps. To confirm this expectation we calculated masses of glue lumps 
associated with static charges in the fundamental representations. In 4-dimensional G2-gluodynamics we found Casimir scaling 
for charges in the representations 7, 14, 27 and 64, similarly as we did in 3 dimensions, although the uncertainties are of course 
larger But within error bars we see no violation of Casimir scaling and this confirms the corresponding results in [16], obtained 
with a variant of the smearing procedure. To see the expected string breaking in 4 dimensions one would need larger lattices 
than those used in the present work. 
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